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Abstract 

The knowledge of the isotropic correlation function of a plane figure 
is useful to determine the correlation function of the cylinders having the 
plane figure as right-section and a given height as well as to analyze the 
out of plane intensity collected in grazing incidence small-angle scattering 
from a film formed by a particulate collection of these cylinders. The 
correlation function of plane polygons can always be determined in closed 
algebraic form. Here we report its analytic expression for the case of a 
triangle. The expressions take four different forms that depend on the 
relative order among the sides and the heights of the triangle. 
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1 Introduction 

In small-angle scattering (SAS) theory samples are conceived of being made 
up of two/three homogeneous phases [IJ In this way the collected scatter- 
ing intensities are used to determine few structural parameters of the observed 
samples. In principle the procedure is very simple. Restricting ourselves for 
definiteness to the case of a two-phase particulate sample and assuming that 
the particle shape is unique, then the effective scattering density of the sample 
reads 

n(r) = (m-najJ^ft-^r/DjO-B,-), (l) 
4=1 

where (rii — n?) is the contrast and the sum is performed over the N particles 
of the samples. Moreover, Dj, R 3 and Clj respectively denote the size, the 
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mass-centre position of the jth particle and the three Euler angles specifying the 
rotation which the jth particle has undergone to. IZ^ denotes the corresponding 
rotation matrix. Finally, p(r) is the function that specifies the particle shape 
since it is equal to one or zero depending on whether the tip of r falls inside 
or outside the particle of unit size (D = 1) with its centre (of mass) set at the 
origin. The observed scattering intensity I(q) [with q = |q| = (4tt/\) sin(6'/2), 
A denoting the ingoing beam particle wave-length and 9 the scattering angle] 
is simply proportional to |n(q)| 2 , the square modulus of the Fourier transform 
(FT) of Eq. (TTJ). The fit of the observed I(q) to |n(q)| 2 , calculated by the 
procedure sketched below, allows one to fix the values of the parameters present 
in Eq. (TTJ). In practice, before proceeding to the best fit analysis, it is necessary 
to exploit any information on the physical structure of the sample. To this aim 
one first considers the convolution of n(r) and then uses the property that its 
FT yields the scattering intensity. The convolution involves the double sum 
Y^,j i—%Tj * Ti where Tj denotes the jth addend of Eq. (TTJ). In handling the 
sum, it is convenient to separate isotropic from anisotropic samples. In the 
first case, the sum of the terms involving particles with different orientations is 
assumed to average to zero provided we substitute the remaining convolutions 
with their isotropic components. These, in turn, are obtained by suitably scaling 
the expression 



which defines the isotropic (auto) correlation function (CF) of the particle of 
unit size and volume V. The FT of 73 (r) is the particle (isotropic) form- factor. 
The scattering intensity is 



I(q) cx ]T D*UD 3 q) + £ D|A 3 73(^-9)73(Ag) e - iq (R ^ Ri) , (3) 



where the tilde and the over-bar respectively denote the FT and the complex 
conjugate. Further approximations as the decoupling pQ and the local monodis- 
perse one [3J allow us to simplify somewhat more expression ([3]). In any case, 
the final expression requires the knowledge of 73(^(7) which is algebraically 
known only for the spherical shape. It can numerically be evaluated by a one- 
dimensional integral for the parallelepiped [H [5], the circular cylinder [5], the 
rotational ellipsoid[7] and the tetrahedron [8] because the CFs of these shapes 
are known in terms of transcendental functions. In all the other cases it requires 
a 3-5 dimensional integral, a condition that makes any numerical analysis with 
other particle shapes much more awkward. Recently it has been shown that 
the CF of a cylinder of arbitrary right section is related to the CF of the plane 
section by a generalized Abel transform [21 HH]- Thus, if one algebraically knows 
the correlation function of a plane figure, the CF of the associated cylinder is ob- 
tained by a quadrature and its FT by a further integration. Generalizing the fact 
that the CFs of the regular triangle, square[T2"]. pentagon[T5] and hexagon [Tl] 
are algebraically knowiiU, Ciccariello[TT] has recently shown that the CF of any 
plane polygon can be determined in closed algebraic form. By this result it is 

x Very recently it has been determined the algebraic expression of the CF of a regular 
polygon whatever the side number 15 . 
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now practically possible to consider cylindrical shapes more involved than the 
circular ones, see e.g. |16) . 

Surprisingly the algebraic expression of the CF of the simplest polygonal shape, 
namely that of a generic triangle, as yet is not known. In this paper we fill 
up the gap. This will be done in the next section. Before concluding this sec- 
tion we wish to emphasize that the mentioned property of the CF of a plane 
polygon can also be used in the case of some anisotropic samples. In fact, in 
the case of fiber samples (see, e.g., Ref.|17j). one uses the anisotropic CF of a 
circular cylinder. By the aforesaid result it becomes now possible to consider 
the correlation function of a polygonal cylinder angularly averaged over the only 
directions perpendicular to the cylinder axis. This result can also be applied 
to the out of plane diffuse component of the intensity collected in grazing inci- 
dence small-angle scattering (GISAS) experiments [T51 |2T)] when the specular 
region is avoided so as to make the kinematical approximation satisfactory. See, 
in particular, the discussion reported in sect. III. A of [T§] . 

Before passing to the illustration of the results relevant to the triangle, we 
recall three basic relations worked out in Ref.[TT] to which one should defer 
for details. They generalize those of the three dimensional(3D) case and are 
respectively related to the values at the origin of the first, second and third 
derivative of the 2D CF. They read 

j'(0) = -L/nS, (4) 
1 M 

7 (0 + ) = ^^[(7r-^)cot^ +1], (5) 

i— 1 

and 

Here L and S are the perimeter length and the surface area, the ipiS are the 
angles at the possible M vertices of the boundary of the figure and R curv (r) is 
the curvature radius at the point of the boundary with curvilinear coordinate 
r. 



2 The triangle case 

According to Ref.[TT] the distances that determine the intervals where 7(f) 
takes an appropriate expression are the lengths of the sides and the heights of 
the considered triangle. Hence we must first establish the order existing among 
these distances. To this aim, given a triangle, we convey to name a the shortest 
side, c the longest one and b the third one so as to always have a < b < c. 
Vertices opposite to sides a, b and c are named A, B and C, and the angles at 
vertices A, B and C as a, (3 and 7 (the context will avoid any confusion with 
the CF's symbol 7). The height from A to a is h a and similarly for hi and h c 
(see figure 1). Owing to condition a < b < c the following inequalities 

a < b < c, h c < h b < h a and a < f3 < 7 (7) 

always hold true. To order the six lengths: a, b, c, h ai hf, and h c we need to know 
if a < h a or a > h a and if the triangle is acute or obtuse. Hence the four distinct 
shapes shown in Fig. 1. 



3 




Figure 1: The triangle's four distinct cases A, B, C and V. Red continuous sides 
— > shape A, green dashed sides — > shape B, blue short-dashed sides — > shape 
C and magenta dot-dashed sides — > shape T>. The side lengths of the shown 
A, B, C and V triangles are (1,1.5,1.611), (1,1.06,1.127), (1,1.5,2.239) and 
(1, 1.5, 2.470), respectively. 



Their names respectively are: 

A -t h c < h b < a < h a < b < c, 7 < tt/2 (8) 

B -> h c < h b < h a < a < b < c, 7 < tt/2 

C — > h c < hi, < a < h a < b < c, 7> tt/2 

V h c < h b < h a < a < b < c, 7> 7r/2. 

These inequalities show that we have six r-intervals. In case A they are: 
[0, h c ], [h Cl hb], [hb, a], [a, h a ], [h a , b] and [b, c]. They will be referred to 
as I A, II A, III A, IV A, V A and VIA, respectively. The same convention is 
adopted in cases B, C and V. Since intervals I A, IB, IC and IV coincide we 
shall simply speak of interval I. The same happens for intervals II and VI. In 
each of the subcases listed in Eq.© we can apply the procedure expounded in 
ref. [TT] and determine the expression of the second derivative of the (isotropic) 
CF of the considered triangle. The procedure is straightforward though some- 
what cumbersome. After putting 

A(a, (3, 7) = 3 + (n — a) cot a + (ir — (3) cot (3 + (n — 7) cot 7, (9) 

B(a,/3,j) = (7r-2a)cotQ; + (7r-2 / 3)cot/3 + 2(l+7Cot7), (10) 

C(a,(3,j) =4 + 2(7r-2a)cota + 7rcot/3 + 7rcot7, (11) 

c(a, /3) = cot a + cot (3 (12) 

and 

^2(2;) = —2xyl — x 2 — 2arccos(a;), (13) 
the final expressions of 7 (r) are 

l" T {r) = A{a,l3, 1 )/{2T:S), < r < h c ; (14) 

7 "„(r) = [c(a,f3)n 2 (h c /r)+A(a,l3, 1 )]/(27rS), h c < r < h b ; (15) 
» c(q, P)n 2 (h c /r) + c(a, j)n 2 (h b /r) + A(a, (3, 7) 
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c(a,(3)n 2 (h c /r) +c(a,j)n 2 (h b /r) + C{a,f3,j) 
4^5 



, a < r < h(fX7) 



[c(a, p)(h(h c /r) + 2c(/3, j)Cl 2 (h a /r) + c(a, j)n 2 (h b /r) (18) 

+C(a,/3, 7 )]/(4^5), h a <r<b- 
l" IIIA {r), h b <r<h a ; (19) 

[c(a, P)n 2 (h c /r) + c(/3, j^iK/r) + c(a, j)n 2 (h b /r) (20) 
+-4(a, 7)] /2wS, h a <r < a; 



a < r < 6; 
h b < r < a: 
a < r < h a ; 
h a < r < 6; 

h b <r < h a ; 
h a < r < a; 
a < r < 6: 



(21) 
(22) 
(23) 
(24) 
(25) 
(26) 
(27) 



c(a, + cQ8, 7 )»2(feg A) + g(a, ft 7) 



, 6 < r < c(28) 



From these expressions the calculation ol 7 (r) is straightforward after noting 
that the primitive of Q 2 (h/x) is function f2i(/i, x) reported below 



Cl±(h, x) 



2/iVx 2 - /i 2 



2xarccos(ft/x). 



(29) 



In fact, integrating Eq. (|28p over the interval [r, c] with 6 < r < c one finds that 

j' VI (r) = [B(a,j3,j)(r-c) + c(a^)(Q 1 (h b ,r)-n 1 (h b ,c)) + 
c08,7)(fii(ft a ,r) - fiiC^.c))]/^, & < r < c, 

where we used the property that 7 (r) is continuous at r = c and therefore equal 
to zero there. It is now convenient to put 

fti(ft„,r, 6) = ni(ft OJ r) - 0i(/i o , b). (30) 

The previous derivative expression then reads 

~f' v A r ) = [B(a,^,j)(r- c) + c(a,-y)Cli(h b ,r,c) + c(^,j) X 

Ai(^ B ,r,c)]/(47r5), 6 < r < c, (31) 

In the same way, integrating Eq. (|18[) over the interval [r, b] with h a < r < b 
and imposing the continuity at r = 6, one finds that 



7 v A r ) 



1 



4tt5 



[2c(/3, 7 )fii(/i a ,r, 6) + c(a,7)fi 1 (/i 6 ,r, 6) + 



(32) 



3(a,0)fti(/i c ,r,&)+C(a, j 0,7)(r-&)] +V V x( & )> h a < r < b, 
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Figure 2: (a): plots of j(r) for triangle A (continuous red curve), B (green 
dashed curve), C (blue short-dashed) and T> (magenta dot-dashed). The con- 
vention in the curve drawing is the same of Fig. 1; (b): plots of 7 (r) (note, 
however, that 7 (r) is scaled by a factor 1/5 in case T>); (c): plots of 7 (r) (the 
shown 7 (r)s are scaled by factors 1/5 and 1/15 in cases C and T>, respectively); 
(d): plots of 7 (r) (the shown 7 (r)s are scaled by factors 1/20, 1/5, 1/10 and 
1/20 in cases A, B, C and D). 



4tt5 



and then, step by step, that 

7' ' 1A ir) 



[c(a,i)Sl\{hb,r,h a ) + c(a, /3)fli(h c ,r,h a ) + (33) 
C(a,/3,7)(r - h a )] +7 VA (h a ), a < r < h a ; 
—g[c(a,j)£li(h b ,r,a) + c(a, (3)Qi(h c , r, a) + 
A(a,/3,j)(r-a)]+~/' IVA (a), h b < r < a; (34) 



:(a, (fe c , r, 6) + -4(a, /3, 7)(r - h b ) 



2nS 



(35) 



' III A 

A(a,/3,j)(r-h c ) 
2nS 



h c < r < h b ; 
+ 7„.(M. 0<r</j c . (36) 
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The expressions of 7 (r) are similarly obtained in the other three cases and are 
reported in the appendix. In the same way, from the resulting expressions of 
the first derivative, using the fact that a primitive of (fii(h, x) is 



Vt{h : x) = — nx /2 + 3h\/x 2 ~h 2 + (x 2 + 2h 2 ) arcsin(/i/a;), 



(37) 



starting from the outermost r-interval and imposing the continuity at the out- 
ermost end points one obtains the algebraic expression of the CF of a generic 
triangle. For instance, the CF expression in the outermost interval is 

7v/ (r) = [s(a,/3, 7 )(r~c) 2 /2 + C (a, 7 )[0(/ lb ,r)-0(/ lb , C )-0 1 (^,c) x 

(r - c)} + c(/3, 7)[fi(ft„, r) - fl(h a , c) - n x (h a , c)(r - c))]] /4ttS. (38) 
This expression, after putting 

tl{h a , r, b) = Q(h a , r) - Q(h a , b) - fl^ha, b)(r - b), 

simplifies into 

7vj (r) = [6(a,/3, 7 )(r-c) 2 /2 + c(a,j)Cl(h b ,r,c) + 



(39) 
(40) 



c(f3 7l Mh a ,r,c)]/(4nS), 
In the remaining interval, for case A one finds 



b < r < c. 



l IV A r ) = 
ll,I A { r ) = 



1 



4tt5 



C{a^ n ){r-bf 



c(a, P)fl(h e , r, b) + c(a, iMh, r, 6) + (41) 



2c(/3, 7 )n(ft a , r, b) + 7 VI (b) (r - 6) + 7v, (&) - 



h a <r < b, 



1 



C(a,/3 )7 )(r-/i ) 



Q(/i CJ r, /i a ) 



2 



c(a, 7)n(/i6, r, feo) + c(a, /3) x 
(ft, a )(r - /i a ) + 7 t , A (/i a ), a<r<h ai 



(42) 



1 



4tt5 



.4(a,/3, 7 )(r-a) 2 



c(a, 7)f2(/i b , r, a) + c(a, f3)£l(h c , r, a) (43) 



+7 IVA (a)(r-a)+'y IVA (a), 



hh < r < a, 



2tt5 



C(a,/3, 7 )(r-^) 2 



+ c(a,^)A(ft c ,r, h b ) 



(44) 



7 UIA ( h b)(r - h b ) +J IIIA (h b ) 

^l(a,/3,7)(r-fe c ) 2 
4tt£ 



h c <r < hb, 

+ 1 IIA { h c){r-h c ) + llIA {h c ), 0<r<h c . (45) 



The expressions relevant to the other cases are given in the appendix. Figures 
2a-d show the plots of 7(f) and its first three derivatives for the four triangle 
cases. The CFs look each other rather similar because the different supports 
cannot be appreciated on the figure scale. The first derivatives show greater 
difference since their behaviours about r = signal that the specific length (i.e. 
the ratio perimeter/surface, see Eq.Q) increases as one passes from case A to V 
(note that the magenta curve is scaled). The differences become more evident 
for the second derivatives and dramatic for the third one. Here the cuspids 
signal algebraic singularities cx (r — /i) -1 / 2 and they involve the only heights 
which meet the opposite side within the triangle. This explains why Fig. 3d 
shows three spikes in cases A and B and only one spike in cases C and T>. 
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3 Conclusions 



We have briefly described the cases where the algebraic knowledge of the CF of 
a plane figure can usefully be applied to analyze observed scattering intensities. 
In particular we have reported here the explicit expression of the CF of a general 
triangle. Besides their intrinsic interest for the shown singularities, which are 
important from the point of view of the figure reconstruction, these expressions 
allow us to evaluate the CF of the associated cylinder by a simple quadrature. 



A Expressions for cases B, C and V 

For completeness we give the expressions of the CF and its fst derivative in 
remaining cases B, C and V. For the first derivative one finds 



B, 



7 V /( r ). 
f 



2nS 



b < r < c; 

2c{/3, 7 )Oi {h a ,r,b) + c{a, 7 )f2i {h b , r, b) + c{a, (3) x 

a < r < b; 

c{a,fi)&i{h c ,r,a) + c{/3, 7)^1 {h a ,r, a) + c{a,j) 



AirS 

fM/ic r, b) + {r- b)C{a, p, 7)] + j' VI {b), 
1 



tti{h b ,r,a) +A{a,P,j){r - a) +"f VB {a), h a < r < a; 



c(a, / 8)Oi(/i c ,r, h a ) + c{a,j)Cli{h b ,r, h a ) + 



2nS 



A{a,f3,j){r - h a ) + J IVB {h a ), h b < r < h a ; 

c{a, f3) (fti(/i c , r) - fti(/i c , h b )) + A{a, (3, j){r - h h ) 



2irS 



A{a,P,j){r-h c ) 
2irS 



7 V7 0), 
f 

4^5 



h c <r < h b ; 
< r < h c ; 

b < r < c; 

[c(a, 7 )fii(/i 6 , r, b) + c{a, jS)fii(/i c , r, b) + C{a, (3, 7 )(r - 6)] 



+7 V/C (6), 



h a <r <b; 



■^-g[c{a,j)Qi{h b ,r, h a ) + c{a, /3)£li{h c ,r, h a ) + 
C{a,/3,j){r-h a )]+j' vc {h a ), a < r < h a ; 



c(a,/3)fSi(fe c ,r,a) + A{a, /3,-y){r - a) 
2ttS 

c(a,/3)Oi(/i c ,r, hb) + A{a, (3,-f){r - h b ) 



(46) 
(47) 

(48) 

(49) 

(50) 
(51) 



A{a,p,-y){r-h c ) 
2nS 



2irS 

+ nf' Iie {h c ), 



+ 7 7V , c ( a ), h b <r<a; 
+ 7 J7C (M, h c <r<h b ; 
< r < h c ; 



(52) 
(53) 

(54) 

(55) 
(56) 
(57) 
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in case C, and 
7v IO W = 7' (r), 



7 VD M = 



7™W 



1 



6 < r < c; 



[c(a, 7)^1 (/ib, r, 6) + c(a, /3)fii (ft c , r, &) + 



4tt5 

C(a,/3, 7 )(r-&)] + 7^», 
A(a, P, 7)(r - a) + cja, P)&i (h c , r, a) 
2nS 



(58) 
(59) 



7,„x.( r ) = — + 



27rS* 



a < r < b; 

+ 7v (°)' ^< r <«; (60) 

(61) 



7„ D W = 
7 ro W = 



7, VB (M, 

c(a,/3)f2i(/i c ,r, fy,) + .4(a, (3, 7 )(r - h b ) 



A{a,p,j){r-h c ) 



2irS 



2irS 



h b <r < h a ; 
+ 7^(M, h c <r<h b ; (62) 
< r < ft c (63) 



in case V. 

For the CF one finds 



IviB M = 

7vbW = 

= 

7««W = 
7„ B (r) = 

7« M = 
in case £5, 

Ivic ( r ) = 
Ivc ( r ) = 

livcif) = 

line W = 



7wM, 

1 



4ttS* 



6 < r < c 

2c(p, 7)A(/i , r, 6) + c(a, y)£l(h b , r, 6) + c(a, /9)ft(fc C) r, 6) 



+C(a,/3, 7 ) 



(r - &) 2 



+ 7 VI (&)(^-&)+7w( & ), a<r<6; 



1 



c(fi, j)Cl(h a , r, a) + c(a, j)Cl(h b , r, a) + c(a, P)Cl(h c , r, a) 



2irS 

+A(a,p,-y) 
1 



(r - a) 2 



2nS 



+ J VB ( a )( r - a )+JvB( a )i K<r<a; 

{r-Kf- 



c(a,P)Q(h c ,r, h a ) + c(a,^)Q(h b ,r, h a ) + _4(a,/J, 7 )- 



+l IVB ( h a)(r - K) + ~f IVB (h a ), 
c(a, P)£l(h c , r, h b ) + A(a, P, 7) 



(r - h b f 



2ttS 

+l' IIIB ( h b)(r - h b ) + j IIIB (h b ), 
A(a,P^)(r-h c ) 2 



h b <r < h a ; 



h c <r < h b ; 



AttS 



+ l 11K ( h c){r~h c )+ llIB (h c ), 0<r<h c 



1 



4tt5 



b < r < c 

c(a, p)Q(h c , r, b) + c(a, j)Cl(h b , r, b) + C(a, p, 7 )(r - fe) 2 /2 



(64) 
(65) 

(66) 

(67) 
(68) 
(69) 

(70) 
(71) 



+l'(b)(r-b)+ lvl (b), 



h a <r <b; 



c(a, P)Cl(h c , r, h a ) + c(a, 7 )Cl(h b , r, h a ) +C(a, p, 7 )(r - h a ) 2 /2 (72) 



2ttS 



+l vc ( h a){r - K) + 7vc {K), 

c(a, p)Q(h c ,r,a) + A(a, p,-y) 



a < r < h„ 



(r - a) 2 



(73) 



2tt5 



+% vc (a)(r -a)+j IVC (a), 
c(a, (3)Cl(h c , r, h b ) + A(a, (3, 7) 
+l IIIC (h b )(r - h b ) + 

"fine 



(r - h b ) 



2-, 



h b < r < a; 



h c <r < h b ; 



lAr) = AM 'll { g +l IIC {hc){r-h c )+ llIC {h c ), 0<r<h c 



in case C, and 
Ivi-d ( r ) = 

l, V T,{ r ) 
in case V. 



1 



AttS 



c(a, P)Ct(h c , r, b) + c(a, i)Q,{h b , r, b) +C{a, p, 7) 



b < r < c: 
(r - b) 2 



2ttS 



+l vl (b)(r-b)+ 7vl (b), 
c{a, /3)h(h c ,r,a) + A(a, (3,^) 



(r - a) 2 



2ttS 



c(a, P)£l(h c , r, h a ) + A(a, p, 7) 
+l' IVV ( h a)(r - K) +7 IV7 ,{h a ), 
c(a, P)Cl(h c , r, h b ) + A{a, P, 7) 



(r - Kf 



a < r < b: 



h a <r < a; 



h b <r <h a 



(r - h b f 



2ttS 

+% 1I Jh b )(r-h b ) + 

IlII'D 

A{a,M)(r-h c ) 2 
AttS 



h c <r < h b : 
+ % rD (h c )(r -h c )+-y Ilr ,(h c ), < r < h c 
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